Trajectories and Stability Regions of the Lagrangian 
Point Li in the Generalized Chermnykh-Like Problem 



Badam Singh Kushvah 

Department of Applied Mathematics, Indian School of Mines 
Dhanbad-826009, INDIA, Phone:+91-326-2235765, Fax:+91-326-2296563 



Abstract 

The Lagrange point Li for the Sun-Earth system is considered due to its spe- 
cial importance for the scientific community for the design of space missions. 
The location of the Lagrangian points with the trajectories and stability re- 
gions of Li are computed numerically for the initial conditions very close 
to the point. The influence of belt, effect of radiation pressure due to Sun 
and oblateness effect of second primary(finite body Earth) is presented for 
various values of parameters. The collinear point Li is asymptotically stable 
within a specific interval of time t correspond to the values of parameters 
and initial conditions. 
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1. Introduction 

The circular restricted three body problem is modification of the three 
body problem where the third body is assumed to have very small mass which 
is infinitesimal in comparison to other two finite masses are called primaries. 
The restricted three body problem is generalized to include radiation pres- 
sure, oblateness of the second primary and influence of the belt. Further the 
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primary bodies are moving in circular orbits about their center of mass. The 
well-known five equilibrium points(Lagrangian points) that appear in the 
planar restricted three-body problem are very important for astronautical 
applications. The coUinear points are unstable and the triangular points are 
stable Szebehely jlj. In the Sun- Jupiter system several thousand asteroids, 
collectively referred to as Trojan asteroids, are in orbits of triangular equilib- 
rium points. But collinear equilibrium points are also made linearly stable 
by continuous corrections of their orbits("halo orbits"). In other words the 
collinear equilibrium points are metastable points in the sense that, like a 
ball sitting on top of a hill. However, in practice these Lagrange points have 
proven to be very useful indeed since a spacecraft can be made to execute a 
small orbit about one of these Lagrange points with a very small expenditure 
of energy Farquhar . Because of the its unobstructed view of the Sun, 
the Sun-Earth Li is a good place to put instruments for doing solar science. 
NASA's Genesis Discovery Mission has been there, designed completely us- 
ing invariant manifolds and other tools form dynamical systems theory. In 
1972, the International Sun-Earth Explorer (ISEE) was established , joint 
project of NASA and the European Space Agency(ESA). The ISEE-3 was 
launched into a halo orbit around the Sun-Earth Li point in 1978, allowing 
it to collect data on solar wind conditions upstream from the Earth Farquhar 
et al. jif. In the mid-1980s the Solar and Heliospheric Observatory (SOHO) 
Domingo et al. jsj is places in a halo orbit around the Sun-Earth Li position, 
about a million miles the Sun ward from the Earth. They have provided 
useful places to "park" a spacecraft for observations. 

The Chermnykh's problem is a new kind of restricted three body problem 
which was first time studied by Chermnykh j6| . This problem generalizes two 
classical problems of Celestial mechanics: the two fixed center problem and 
the restricted three body problem. This gives wide perspectives for applica- 
tions of the problem in celestial mechanics and astronomy. The importance 
of the problem in astronomy has been addressed by Jiang and Yeh [7| . Some 
planetary systems are claimed to have discs of dust and they are regarded 
to be young analogues of the Kuiper Belt in our Solar System. If these 
discs are massive enough, they should play important roles in the origin of 
planets'orbital elements. Since the belt of planetesimal often exists within 
a planetary system and provides the possible mechanism of orbital circular- 
ization, it is important to understand the solutions of dynamical systems 
with the planet-belt interaction. Chermnykh's problem has been studie d by 
many scientists such as Papadakis [sl], Jiang and Yeh jif, Yeh and Jiang 10 
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Papadakis and Kanavos and reference their in. 

The goal of present paper is to investigate the nature of colhnear equihb- 
rium point Li because of the interested point to the mission design. Ahhough 
there are two new equihbrium points due to mass of the beh(larger than 0.15) 



Jiang and Yeh [12|, Yeh and Jiang [101 but they are left to examine. All the 
results are computed numerically with the help of computer because pure 
analytical methods are not suitable. The actual trajectories and the stabil- 
ity regions of Li however is more complicated than the discussed here. But 
for specific the time intervals, and initial values, these results provide new 
information on the behavior of trajectories around the Lagrangian point Li 
for different possible set values of the parameters. 

2. Location of Lagrangian Points 

It is supposed that the motion of an infinitesimal mass particle is influ- 
enced by the gravitational force from primaries and a belt of mass Mi,. The 
units of the mass, the distance and the time are taken such that sum of 
the masses and the distance between primaries are unities, the unit of the 
time i.e. the time period of mi about 777-2 consists of 27r units such that the 
Gaussian constant of gravitational = 1. Then perturbed mean motion 
n of the primaries is given by r?,^ = 1 + + ^ 2^,^2^3/2 ; where T = a + b, 
a, b are flatness and core parameters respectively which determine the den- 
sity profile of the belt, = (1 — fJ.)ql^^ + ^2 = is the oblateness 
coefficient of 7772; r^, Tp are the equatorial and polar radii of m2 respectively, 
r = + 7/^ is the distance between primaries and x = fi{t),y = f2{t) are 
the functions of the time t i.e. t is only independent variable. The mass 
parameter is u = ^5^(9.537 x 10"^ for the Sun- Jupiter and 3.00348 x 10"^ 

" mi+m2 ^ 

for the Sun-Earth mass distributions respectively ), gi = 1 — ^ is a mass 
reduction factor and Fp is the solar radiation pressure force which is exactly 
apposite to the gravitational attraction force Fg. The coordinates of t?7i, 
7772 are (— /7,0), (1 — ^,0) respectively. In the above mentioned reference 



system and Miyamoto and Nagai [13|| model, the equations of motion of the 



infinitesimal mass particle in the xy-pla.ne formulated as [please see Kushvah 
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X — 2ny = Q^, (1) 
y + 2nx = ^y, (2) 



3 



where 



2 (1 - iJ,)qi(x + ^) fi(x + ^-l] 
n X ^ ^ 



3 iiA2{x + fi — 1) 
2 r| 



MbX 



(r2 + r2) 



3/2 



2 /U?/ 3/1^2?/ 

n y - ■ 



r 



2 ri 



(r2 + T2) 



3/2 



2 ri r2 2r| (^2 + 7^2)1/2 

ri = + /i)2 + y"^, r2 = ^/ {x + - 1)^ + y'^. 

From equations ([1]) and ([2]), the Jacobian integral is given by: 

E =^(^x'^ + y^^ - Q{x, y, x, y) = (Constant) (4) 

which is related to the Jacobian constant C = —2E. The location of three 
collinear equilibrium points and two triangular equilibrium points is com- 
puted by dividing the orbital plane into three parts Li,L4(5): fi < x < 
(1 — /i), L2: (1 — //) < X and L3: X < —fi. For the collinear points, an 
algebraic equation of the fifth degree is solved numerically with initial ap- 
proximations to the Taylor-series as: 

^(L,) = + + (5) 

xiU) = l + + + (6) 

^ ^ 12 20736 248832 ^ ^ 

(8) 

The solution of differential equations ([1]) and ([2]) is presented as interpolation 
function which is plotted for various integration intervals by substituting 
specific values of the time t and initial conditions i.e. x(0) = x(Lj),?/(0) = 



4 




Figure 1: The position of equilibrium points when T=0.01, qi — 0.75, A2 = 0.05 and 
Mb — 0.4, panel (I):Red doted curves and blue points for Sun- Jupiter mass distribution, 
blue curves and black points for Sun-Earth mass distribution, (II): Position of Li, L2 with 
respect to Jupiter's and Earth's is shown in zoom 

where i = 1 3 and a;(0) = | — /i, ?/(0) = ±^ for the triangular equihbrium 

points. The equihbrium points are shown in figure [1] in which two panels 
i.e. (I) red solid curves and blue points correspond to the Sun- Jupiter mass 
distribution and blue dashed curves and black points correspond to the Sun- 
Earth mass distribution. Panel (II) show the zoom of the neighborhood of 
Li,L2. The numerical values of these points are presented in Tabled] One 
can see that the positions of Li , L3 appeared rightward and the positions of 
L2,L4 ( L4 is shifted downward also) are shifted leftward in the Sun-Earth 
system with respect to the position in the Sun- Jupiter system. The nature 
of the L5 is similar to the L4. The detail behavior of the Li with stability 
regions is discussed in sections. [3]&|11 



3. Trajectory of Li 

The equations (II]l2]) with initial conditions x{0) = x{Li), y{0) = 0, x'{0) = 
y'{0) = are used to determine the trajectories of Li for different pos- 
sible cases. The origin of coordinate axes is supposed to the equilibrium 
point at time t = to draw the figures which show the trajectories of 



5 



Table 1: Location of equilibrium points when T=0.01, qi = .75, A2 = .05 and Mi, = 0.4 



Sun- Jupiter 

Li X y 


Sun-Earth 

X y 


Li 0.774577 
L2 1.09493 
L3 -0.786195 
L4 0.410603 0.669308 


0.78569 
1.0232 
-0.785732 
0.393072 0.680342 



the point in consideration. They are shown in figure [2] with six panels 
i.e the panels (I-IIl) show the trajectory moves about the origin (Li at 
t = 0) with X e (0.990093,1.00916) , y E (-0.0061448,0.00587171), the 
energy E G (-12706.5(t = 22.66), -5.08226(t = 0)) and the distance r{t) e 
(0.990093(t = 0), 1.00916(t = 55)). The panels (1-111:127 < t < 129.6) show 
the trajectory moves away from the origin (Li at t = 0) after a certain value 
of the time t, with x E (0.990093,1.00916) , y E (-0.0061448,0.00587171), 
minimum energy E = —1447 found at the time t = 128.52, and the energy 
> for t > 128.88. 

Figure [3] is plotted for qi = I, Mb = and A2 = 0.05 with six panels 
(1-111:0 < t < 0.06) and (1V-V1:0.06 < t < 1) which describe the effect 
of oblateness of Earth to the trajectory of Li. The graphs plotted against 
time which describe behavior of trajectories to equilibrium points not the 
point itself is moving with time, x = —1.91954 x 10~^^(f = 0.06) to x = 
0. 99405 (t = .05) coordinate y is deceasing function of time that reach maxima 
-0.0000530614, at t = 0.04 and minima -0.000105662 at time t = 0.05 again 
it deceases and reach at value —2.5677 x 10^''(t = 0.06) . Initially energy has 
negative values for time << t < 0.059 decreases with time t which attains 
minimum value -2.64032 x 10^{t = 0.059) then strictly increasing function 
that attains positive values after time t = 0.0594. In time interval (0.2, 
0.6) energy one time returns down that again it tend to very large (infinite) 
positive value. It is clear from panels (IV- VI) the trajectory move far from 
the Lagrangian point Li after time t = 0.0594. The distance r{t) from this 
point to the trajectory is increasingly periodic for time < t < 0.6 then tend 
to very large. 

The effect of radiation pressure, oblateness and mass of the belt is con- 
sidered in figure 4, panels (I&III) describe the trajectory and panels (II&IV) 
show the energy with respect to the time t. The mass reduction factor 



6 




Figure 2: The Panels (I-III):0 < t < 128.23 and (IV-VI):127 < t < 129.6 in which (I and 
II) show the trajectory of Li, (II and V) show energy- versus time and (III- VI) show the 
local distance of trajectory at time t form the initial points i.e. t = the other parameters 
are T=0.01, = 1, ^2 = and Mi, = 0. 
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Figure 3: The Panels (I-III):0 < t < 0.06 and (IV-VI):.06 < t < 1 in which (I and II) 
show the trajectory of Li, (II and V) show energy- versus time and (III- VI) show the local 
distance of trajectory at time t form the initial points i.e. t = the other parameters are 
T=0.01, qi = 1, A2 = 0.05 and Ah = 0. 
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Table 2: Trajectory of Li and the energy when T — 0.01, qi = .75, Mb = 0.2 



A2 


time t 


X 


y 


Energy E 


0.25 


0.000 


0.990093 


-1.27593 X 10^^^ 


-3946.49 




0.002 


0.990333 


-4.4421 X 10^^ 


-4460.93 




0.004 


0.991106 


-3.63139 X 10-6 


-6764.39 




O.OOd 


0.99/646 


— 1.2682^ X iO 


-17501.7 




0.008 


0.996301 


-285308 X 10-5 


-544626. 




0.010 


4.03799 X 10^^ 


-1.06467 X 10^^ 


7.84745 X 10^1^ 


0.50 


0.000 


0.990093 


7.51113 X 10-33 


-7887.36 




0.002 


0.99058 


-9.7398 X 10-^ 


-10146.9 




0.004 


0.992298 


-8.12979 X 10-*^ 


-27763.1 




0.006 


-8.35666 X lO^s 


-1.70025 X 10^^ 


2.17394 X 10^°^ 




0.008 


-3.05605 X 10^6 


-6.21786 X 10^^ 


2.74317 X 10119 




0.010 


-9.05111 X 10^7 


-1.84154 X 10^*^ 


6.20464 X 10121 


0.75 


0.000 


0.990093 


-1.80556 X 10-34 


-11828.2 




0.002 


0.990836 


-1.58382 X 10-6 


-17463.6 




0.004 


0.993821 


-1.33765 X 10-5 


-125336. 




0.006 


6.67303 X 10^5 


1.19389 X 10^4 


4.20091 X lOiis 




0.008 


1.02535 X 10^® 


1.83448 X 10^6 


1.32379 X 10122 




0.010 


1.19678 X 10^9 


2.1412 X 10^7 


6.74915 X 10123 



q\ = 0.75 and Mh = 0.2 are taken to plot the graphs variation in values of 
these parameters have similar effect. In the panels, solid blue lines represent 
A2 = 0.25, red dashed hues correspond to A2 = 0.50 and doted black hues 
for A2 = 0.75. One cane see that the trajectory move very far from the Li 
the energy is positive after a certain value of the time t. Details of trajectory 
and energy is presented in Tabled for various values of parameters. One can 
see that x is an increasing function of the time but y is an initially decreasing 
function for certain values of the time, then it becomes a strictly increasing. 
Similarly the energy E is negative and went downward but after some spe- 
cific the time for each cases it becomes positive and strictly increasing and 
attains very large positive value. 
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Figure 4: The Panels (I-III):0 < t < 0.005 and (II-IV):.00 < i < 1 in which (I and II) 
show the trajectory of Li, (II and V) show energy- versus time and (III- VI) show the local 
distance of trajectory at time t form the initial points i.e. t = the other parameters are 
T=0.01, gi = 1, ^2 = 0.05 and Mb = 
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4. Stability of 



Suppose the coordinates {xi,yi) of Li are initially perturbed by changing 

x(0) = xi + ecos(0),y(O) = Hi + esin(0) where (p = arctan ^ij^pl^) ^ 

(0,27r),0 < e = ^(x(0) - xi)^ + {y{0) - yif < 1. The cj) indicates the 
direction of the initial position vector in the local frame. If the e = means 
there is no perturbation. It is supposed that the e = 0.001 and the = | to 
examine the stability of Li. Figure E] show the path of test particle and its 
energy with four panels i.e. the panels (I&III):gi = 0.75, 0.50,^2 = 0.0, in 
(I) trajectory of perturbed Li moves in chaotic-circular path around initial 
position without deviating far from it, then steadily move out of the region. 
In (111) the test particle move in stability region and returns repeatedly on 
its initial position. The blue solid curves represent = 0.25 and dashed 
curves represent = 0.50. It is clear form panel (III) that bounded region 
for Mb = 0.25 is t < 2500 and for = 0.50, t < 2600. 

The effect of oblateness of the second primary is shown in figure |6] when 
qi = 0.75, Ml, = 0.25. The panel (I) shows the trajectory of perturbed point 
Li and (II) shows the energy of that point. The blue doted lines correspond 
to Ai = 0.25 and red lines for A2 = 0.50. One can see that the oblate effect 
is very powerful on the trajectory and stability of Li. When = 0.0 the Li 
is asymptotically stable for the value of t which lies within a certain interval. 
But if oblate effect of second primary is present (^42 7^ 0), the stability region 
of Li disappears when this effect is increases. Further all the results presented 
in the manuscripts are similar to the results obtained by fld\ , Kushvah 17 ■ 



5. Conclusion 

The numerical computation presented in the manuscript provides remark- 
able results to design trajectories of Lagrangian point Li which helps us to 
make comments on the stability(asymptotically) of the point. We obtained 
the intervals of the time where trajectory continuously moves around the Li, 
does not deviate far from the point but tend to approach (for some cases) it, 
the energy of perturbed point is negative for these intervals, so we conclude 
that the point is asymptotical stable. More over we have seen that after 
the specific time intervals the trajectory of perturbed point depart from the 
neighborhood and goes away from it, in this case the energy also becomes 
positive, so the Lagrangian point Li is unstable. Further the trajectories and 
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Figure 6: Effect of oblateness coefficient A2 on tlie stability of Li panel (I) trajectory 
(II) energy of perturbed point Li in which blue solid curves for A2 = 0.25, red curves for 
A2 = 0.50 

the stability regions are affected by the radiation pressure, the oblateness of 
the second primary and mass of the belt. 
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